Josephson effect in superconducting constrictions with hybrid SF electrodes: pecuHar 
properties determined by the misorientation of magnetizations 



O 

o 



> 
o 



o 

-I— > 



I 



■ 

X 



A. V. Zaitsev 

Institute of Radioengineering and Electronics of the Russian Academy of Sciences, 125009 Moscow, Russia 

' Josephson current in SFcFS junctions with arbitrary transparency of the constriction (c) is investi- 

I gated. The emphasis is done on the analysis of the supercurrent dependencies on the misorientation 

■ angle 9 between the in-plane magnetizations of diffusive ferromagnetic layers (F). ft is found that 
the current-phase relation I{ip) may be radically modified with the 9 variation: the harmonic Ji sin ip 
vanishes for definite value of 9 provided for identical orientation of the magnetizations {9 = 0) the 
junction is in the "tt" state. The Josephson current may exhibit a nonmonotonic dependence on the 
misorientation angle both for realization of "0" and "tt" state at ^ = 0. We also analyze the effect 

' of exchange field induced enhancement of the critical current which may occur in definite range of 

(N : 9. 

r—l. PACS numbers: 73.21.-b, 74.45. +c, 74.50.-|-r 

; 

^ The dc Josephson effect in junctions with ferromagnetic interlayers (F) exhibit exchange field induced remarkable 

?H features jTiamely, the transition from a "0" state to a "tt" state of junctions 0,01, enhancement of the critical Josephson 
current 0,01 and others. These features have been intensively studied in recent years theoretically and experimentally 
^ • (see Reviews 00,0 and Refs. therein). In spite of a lot of works devoted to this problem there are several important 
^ ; [ ones which require additional investigation. The purpose of the present work is the study of one of these problems. 
d . Namely we investigate the Josephson current and peculiar properties of SFcFS junctions (c denotes a constriction) with 

■ arbitrary misorientation angle 9 between the in-plane oriented magnetizations of the ferromagnetic metal (F) layers. 

1 , This investigation is important because the supercurrent is very sensitive to the mutual orientation of magnetizations 
' ^ ' and, besides, because the variation of 9 gives additional opportunity to switch experimentally between the "0" and 

"tt" states. In the considered junctions the constriction is an insulator characterized by the arbitrary transparency or 
short (in comparison with the coherence length) diffusive metal channel. For the F layers the dirty limit, i.e. diffusive 
regime of electron transport is assumed. Note that for the case of ballistic transport through the F layers of quantum 
S-FIF-S contacts with arbitrary misorientation angle 9 between the magnetizations in the F layers the Josephson 
current has been studied by Barash et al. . In the present paper we investigate the angle ^-dependence of different 
peculiar properties of the supercurrent in SFcFS junctions. We show that variation of 9 may result in the transition 
' between the "0" and "tt" states and the disappearance of Ii sin (p harmonic of the supercurrent for certain value of 9. 

■ Effect of exchange field induced enhancement of the critical current and its nonmonotonic dependence on the angle 9 
C^l ■ are also studied. 

' As a model of the constriction, we consider an aperture of a small radius in a thin impenetrable screen dividing 
two different electrodes. The constriction is supposed to include a barrier with the transparency D — D{'d), where 
is the angle between the electronic trajectory and normal to the junction plane. Consider first the case of ballistic 
constriction of small size assuming that its radius is small in comparison with Ip, Wi?/A, dp, where Ip and vp, are the 
mean free path and Fermi velocity in the F layers, di;- is their thickness, A is the energy gap in the superconductors. 
Under these assumptions the current may be expressed through the Matsubara Green's functions of the F layers 
^ . Gi 2 as follows 0[lS|l3: 

a : 

8: /(^) = ^^(j(^„)), (1) 

n—O 

, D G_G+ D [G2,Gil_ 

J = -rTrao ® fg . + = -Trao ® 7-3 L ' ^ . , 

4 {l-D)l + DGl 4 2{2-D)l + D[G2,Gi] + 

where w„ = TTT{2n + 1) is the Matsubara frequency, G± = (62 ± Gi)/2, [G2,Gi]± = G2G1 ± G1G2, i?o = 
Att'^H'^ /{e^p^A) is the Sharvin resistance of the contact, A being the area of the contact, angular brackets denote 
angular-averaging ((...)) — 2 (...) cos -dd cos -d, (7^ and are the Pauli {k = 1,2,3) or unite (fc — 0) matrices in 
the spin and particle- hole spaces, respectively, 1 = tq ® ctq As was noted above we consider the case of dirty thin F 
layers, i.e. their thickness is supposed to satisfy the conditions Ip <^ dp {hDp/TcY^^ where Dp is the diffusive 
coefficient in the F layers, is the critical temperature of the superconductors. For such layers the space variation 
of the Green's functions is negligible therefore they are given by the expression G^ = exp{iipjP3/2)gj exp{—iipjP3/2), 
where pa = 73 (8)(To, (fij is the phase of the order parameter of the jth S-electrode and gj is determined by the following 
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equation (for details of the derivation see, e.g., [l^): 



[gj,ujp3 + ihjS + ebgs]+ = 6 



(2) 



where Sb = fi {Db) vp /^dp ~ Dp /2Rb(JFdF^ Db and Rb are, respectively, the transparency and the resistance per unite 
area of the S-F interface, ap is the conductivity of the F metal, gs is the Green's function of the superconductors 
at the S-F interfaces, hi. 2 are the exchange fields in the F layers, the matrix vector S is defined as J^, 14] S = ct 
T+—(J2^(J2 ® T^, where a — (cti, (72, (T3), t± — {tq ± fs)/2. Note that the value of the exchange fields, |hi^2| = ^ is 
supposed to be small with respect to h/rp and the Fermi energy of the F layers, where Tp = Ip/vp; for considered 
here dirty layers h/rp ^ A. We assume for simplicity that the S-F interfaces are identical and that their transparency 
Db is small, therefore the effect of the F layers on the superconductors is negligible. The only difference between the 
F layers is related with the orientations of the exchange fields, hi 2, which being parallel to the layers, make an angle 
6 with each other. In this paper we confine ourselves to the case of conventional s-wave superconductors. Therefore 
to carry out further calculations it is convenient to represent the current with the use of the transformed Green's 
functions Gj = exp(i(pjP3/2)gj ex]3{—iipjP3/2) where 



9 J = UgjU- 



(3) 



U = (<To ® T+ + (T2 ® T-). Taking into account that particle (hole) projection matrices, obey the relation 

TaTf3 — TQ.Saf3 it is casy to check that U is the unitary Hermitian matrix, U = — [/"*". From (0), ^ one can find 
the representation for gj with the help of the projection matrices in the spin space, P±{nj) — ^((Tq ± Hj't), where the 
unite vectors nj — hj/ \hj\ . Taking into account the relation Pa{n)Pi3{n) = Pa{n)Sa0, we get 



9j 



Q = ± 



1 



1 



5a Pa(nj) = -(5+ -I- g-) ® (To + -(.g+ - g-) n^a 



(4) 



where ga = gah + /aT2, 5q + /a = 1; 



J. 1 Ja J, : 



(5) 



— \/<^a + Cua — + icth + £b55(w), Ap — Ebfsi^)- The form of the Green's function representation Q does 
not depend on the specific expressions for g^; it is valid for spatially- homogeneous orientation of the exchange fields 
in the F layers. Note that being written via the transformed matrices Gj, the expression for the current coincides 
with (1). The advantage of the representation 10} is related with significant simplification of the matrix structure 
of the Green's functions which are given by sum of two terms determined by the direct product of matrices in the 
particle- hole and spin spaces. With the help of we get 



[G2, Gi]± = 2 X! 



'_G2a,Gip]+ (E> [Pjn2),Pp{ni)]± + [G2a,Gif3]- (E) [Pa(n2), P,3(ni)] 



(6) 



where Gja = <?ct'?3 + exp(i(^j-f3)/af2. Taking into account that 



[Pa{n2),Pp{ni)]+ = -[1 -I- a/3n2ni 4- (an2+/3ni)(7], 
[Pa(n2),P/3(ni)]_ = ^[n2xni]a. 



after some calculations we obtain the following expression for the current 

SttT sin (/3 / ^i^) 



E 



e^o ^^\Bo-[A^ + A{^)]s\n^{^) 



(7) 



(8) 



where Lp = ip2 ~ '■Pi, A{{p) = Ao + Ai sin^ 



Aq = 2D 



2 cos2 ( 1 j Re /2 + sin2 r-]\f\^{2 + Dq) 



Ai = -AD^\f\\ Bo 



2 + Dsm' [-]q 



, q^\g\' + \f\'~l. 
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here g = g+{ujn), f = f+{^n)- Eq. © is significantly simplified for the cases of parallel (p) and antiparallel (a) 
magnetizations : 

_ STTTsiny^ ^ / Rc {Df/[1 - Df sin^ (f )]) , p 

^^\(2D\f\'/[2 + D{\g\' + \f\'cosv~l)]),a 



Note that for the case 9 = Eq. @ reduces to the result obtained in Ref . ^3 • In general case the misorientation angle 
dependence of the current is rather nontrivial if the transparency D is not small: it does not reduce to the expression 
(if hi^ha) [13 

(^)cos2 (^0 +/W(^)sin2 . (10) 

Simple form for the 6*— dependence pU|l is valid only in the limit of small transparencies (when I'-P'°-'> (ip) = I'^P'"-^ sin ip); 
for this case it was first obtained in Ref.0]. 

For numerical analysis of the Josephson current we apply the obtained results assuming that superconductors are 
described by the BCS theory therefore gsi^) = w/(w^ + A^)^/^, fs{uj) = A/{u)^ + A^)^/^. The results of numerical 
calculations of I{ip) -dependences for a junction with the transparency D = 0.5 are presented in Fig.l for a set of values 
of 6. It shows that the function I{(p) may strongly deviate from the sine/?— dependence and is strongly modified with 
increasing angle 0, especially in the vicinity of misorientation angle 9 = ^o-n- (temperature Tott) which corresponds 
to the transition between the "0" and "tt" states of junction. Note that the model with the angle- independent 
transparency corresponds to the case of quantum single-mode constriction with the resistance Rq = nh/e^. At 
definite 9 = 9* (temperature T*) which with the accuracy higher than 5 percents coincides with the value 6*0^ {Tq-k) 
the harmonic Ii sinip of the current (/(</j) = Ii sin 1^9 -|- I2 shi{2ip) + ...) vanishes; the sign of Ii is changed with the 
variation of 9 (or temperature) near 9* (T*). Living the detailed analysis of supercurrent temperature dependences 
for a separate presentation we show here the typical phase diagram of the junction (Fig. 2) for not too small exchange 
fields (in comparison with £{,). Another interesting property of the supercurrent is nonmonotonic dependences of the 
critical current Ic = max|/((^)| as a function of the angle 9. Note that a similar behavior of the critical current for 
quantum S-FIF-S junctions with the ballistic transport through the F layers has been analyzed in 0. We find that 
in junctions with diffusive F layers the nonmonotonic dependence /c(^) (Fig. 3) may occur for realization of both "tt" 
and "0" state for identical orientation of the magnetizations (unlike the case studied in Q)- The results of numerical 
calculations of the Ic{h) dependences for a low transparency junction at various angle 9 are presented in Fig. 4. It 
shows that the ratio /c(ft.)//c(0) may exceed unity, i.e., the effect of the exchange field induced enhancement of the 
critical current 0,0] may occur in some range of ft, if 6* exceeds definite value (7r/4 < 6* < tt for small D). 

For the case of SFcFS junctions with a short diffusive constriction the following expression for the current may be 
obtained with the use of the Green's function technique ^3 

I{^)= [' p{D)l{^,D)dD, (11) 
Jo 

where l{ip,D) is determined by ijHJl with Rq = Trh/e^ being the resistance of a single mode constriction, and the 
distribution of transparencies is determined by the function p{D) = {Rq/2Rn)/ D{\ — D)^/^ (Rn is the normal state 
resistance of the constriction) which coincides with the density function found by Dorokhov 'lr\. For the cases of 
parallel and antiparallel magnetizations Eqs.®, (|11|) reduce to 

le^^^sin^^ r Re F^VF (1-2/2 sin^d)), p 

where W{x) = (1 — a:^)"^/^ arctan[(l — x)/{l + x)]^^^. Note that in the case ^ = Ea. lfT^ reduces to the result 
obtained in Ref. for h = Eas. (|12|l coincide and in the limit e?, ^ A reproduce the Kulik-Omelyanchuk formula 
[igj i for a diffusive ScS constriction. Numerical calculations show (Fig. 5) that the above discussed peculiar properties 
manifests itself also in junctions with diffusive constriction. 

In conclusion, we have developed a microscopic theory of supercurrent in SFcFS junctions with the arbitrary 
transparency of the constriction and diffusive electron transport in the F layers for arbitrary angle 9 between the in 
plane exchange fields in the F layers. We showed that the current-phase dependence may be strongly modified with 
the variation of 9. If the junction is in the "tt" state for identical orientation of magnetizations {9 = 0), the variation 
of 9 may result in the transition between the " 0" and " tt" states and disappearance of the Ii sin ip harmonic of the 
supercurrent at definite value of 9. Nonmonotonic dependence of the critical current as a function of 9 may occur for 
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both "tt" and "0" state of junction corresponding to the case of identical orientation of the magnetizations if the 
transparency of the junction is not smaU. 
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FIG. 1: Fig.l. Current-phcise relations at various misorientation angle 9 for a SF-I-FS junction with the transparency of the 
barrier D = 0.5 at h/lS. = 1.5, T/A = 0.15, £6 = A. For < 0ot, = 7r/4 (^ovr < 9 < tt) the junction is in the "tt" state 
("0" state). For 9 in the vicinity of 9otv second harmonic /2sin(2i^) of the current dominates, first harmonics Iisimp vanishes 

at ^ = ^ !^ 9oTr- 




FIG. 2: Fig.2. The {T,9) phase diagram of a SF-I-FS junction with the transparency of the barrier D = 0.5 and h/A = 
1.5, Sb = A. Black region corresponds to the "tt" state of the junction. 
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FIG. 3: Fig. 3. Critical current as a function of the misorientation angle 9 normalized to its value at = for a SF-I-FS 
junction with the transparency of the barrier D = 0.5; £(, = A, T/A = 0.01; these dependences correspond to "0" state of the 
junction. 




FIG. 4: Fig.4. Exchange field dependence of the critical current Ic{h) (normalized to /c(0) = Ico) at various misorientation 
angle 9 for SF-I-FS junction with the transparency of the barrier D = 0.1 at T/A = 0.01, £6 = A. 
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FIG. 5: Fig.5. Current-phase relations at various 9 for a diffusive SF-c-FS junction at ft/A = 1.5, T/Tc = 0.1, = A. For 
< 9oiv = 0.647r (^ovr < 9 < t:) the junction is in the "tt" state ("0" state). For 9 in the vicinity of ^ovr second harmonic 
I2 sin(2i^) of the current dominates, first harmonics /i sin<^ vanishes aX 9 = 9* !V ^ott- 



